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Abstract. The detailed investigation of the single pion production reactions pp — ► pmr + and pp — >• pp7r° at 
the incident proton momentum 1628 MeV/c has been carried out. The data are analyzed in the framework 
of the event-by-event maximum likelihood method together with the pp —¥ ppiv data measured earlier in 
the energy region below 1 GeV. At 1628 MeV/c the largest contributions stem from the 3 P2, 3 Pi, 3 fo, 
1 D2 and 3 F2 initial partial waves. 

PACS. 13.75.Cs Nuclon-nucleon interactions 
scattering 



13.85.Lg Total cross sections - 25.40.Ep Inelastic proton 



1 Introduction 

Nucleon-nucleon interaction is one of the most impor- 
tant processes in nuclear and particle physics and exten- 
sively has been studied over a wide energy range. The sin- 
gle pion production in the NN interactions is the main 
inelastic process at energies below 1 GeV. Despite the 
fact that a lot of experiments have been performed, many 
questions on this process are not yet answered. One of 
them is the question about a contribution of the isoscalar 
(I = 0) channel to the inelastic neutron-proton collisions. 
Since the neutron-proton scattering amplitude contains 
both isoscalar and isovector (/ = 1) parts, a detailed in- 
vestigation of the single pion production in the pp col- 
lisions (isovector contribution only) might give the most 
accurate information about the isovector channel which, 
in combination with the neutron-proton data, would allow 
one to extract correctly the contribution of the isoscalar 
channel. 

Various theoretical models, more or less successful, 
were put forward while the data on the pion production 
in the NN collisions were accumulated. Most of them 
are constructed for the energy region near the produc- 
tion threshold and can not be applied at higher energies. 
For the energy range about 1 GeV the one-pion exchange 
(OPE) model [T] assumes a dominance of the one-pion ex- 
change contribution to the inelastic amplitude. Pole dia- 
gram matrix elements were calculated using a beforehand 
form factor. The form factor function was obtained by 
fitting experimental data, so in fact this was a semiphe- 
nomcnological model. In addition, only the P33 partial 
wave was taken into account in the intermediate ttN chan- 
nel [2J. Nevertheless, this model predicts (up to a normal- 
ization factors) with a reasonable accuracy the differential 



spectra of the pp — > pnir + and pp — > ppn reactions in 
the energy range 600-1300 MeV [1[3]. At the same time 
discrepancies between the measured total cross sections 
for these reactions and the model predictions are sizeable 
enough. 

It should be noted that the experimental data on the 
differential spectra of pp — > ppir reaction near the energy 
of 1 Gev are more scarce than the data for pp — > pnn + 
channel. The KEK data [J] contain the information on to- 
tal cross sections only. Our data on the differential spectra 
of the pp — > ppn reaction at the energies 990 MeV and 
900 MeV were published earlier [5]. It would be impor- 
tant to perform an accurate measurement of the differen- 
tial cross sections in the middle of this energy region. Here 
we present our investigation of the pp — > pnn + and pp — > 
ppn reactions at 940 MeV and determine contributions 
from various partial waves to the single pion production 
processes. 



2 Experiment 

The experiment was performed at PNPI 1 GeV synchrocy- 
clotron. The events were registered by a 35 cm hydrogen 
bubble chamber disposed in the 1.48 T magnetic field. 
The proton beam (after corresponding degrader for the 
momentum 1628 MeV/c) was formed by three bending 
magnets and by eight quadrupole lenses. The incident pro- 
ton momentum value was inspected by the kinematics of 
the elastic scattering events. The accuracy of the incident 
momentum value and momentum spread was about 0.5 
MeV/c and 20 MeV/c (FWHM) correspondingly. A total 
of 4.6 x 10 5 stereoframes were obtained. The frames were 
double scanned to search for events due to an interac- 
tion of the incident beam. The double scanning efficiency 
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was determined to be 99.95%. Approximately 8xl0 3 two- 
prong events were used for subsequent analysis. 

The 2-prong events selected in the fiducial volume of 
the chamber were measured and geometrically reconstruc- 
ted. The reconstructed events were kinematically fitted to 
the following reaction hypotheses: 



p + p^ p + p, (1) 

p + p^p + n + 7r + , (2) 

P + P -> P +P + 7T°, (3) 

p + p^d + n + , (4) 

p+p -)• d+ 7T+ + 7T°. (5) 



The identification of the events was performed on the 
X 2 criteria with confidence level less then 1%. If the event 
had a good \ 2 f° r the elastic kinematic (4C-fit), it was 
considered as elastic one. Stretch functions for the three 
kincmatical variables of a track (the inverse of the mo- 
mentum, the azimuthal and dip angles) were examined 
for elastic scattering events to make sure that their er- 
rors and the error for bubble chamber magnetic field were 
properly given. If there was only one acceptable fit for 
the event it was identified as belonging to this hypothesis 
(with a check of stopping ir + track on the presence of the 
7r — > /it — > e decay). If several inelastic versions revealed 
a good x 2 j w e used visual estimation of the bubble den- 
sity of the track to distinguish between proton (deuteron) 
and pion. For few events, even after repeated measure- 
ments, the fit revealed only one acceptable hypothesis but 
with a large \ 2 ■ If such event had x 2 l ess than 50 for the 
4-constraint fit or less than 20 for the 1-constraint fit, it 
was taken into account for the calculation of the cross sec- 
tion value of the process corresponding to this hypothesis 
but not included in differential spectra. 

There were events which failed to fit any hypothesis. 
These no-fit events were investigated on the scanning ta- 
ble, and most of them appeared to be events with a sec- 
ondary track undergoing one more scattering near a pri- 
mary vertex. 

There were also events unfit for the measurements, 
e.g. events with a bad vertex or superimposed tracks. The 
number of such events was counted approximately to be 
7%. The total number of 2-prong events which did not 
pass the measurement and fitting procedures was counted 
to be less than 10%. These unidentified events were ap- 
portioned to the fraction of the fitted hypotheses for the 
total cross section measurements. 

Missing mass distributions for accepted events show 
clear peaks - at zero for the elastic scattering, at the 7r° 
mass squared for single neutral pion production and at 
the neutron mass squared for the process with neutron in 
the final state. 

The standard bubble chamber procedure was used to 
obtain absolute cross sections [3] for the elastic and sin- 
gle pion production reactions. These values together with 
statistic of the events are listed in Tabic [TJ 



Table 1. Numbers of events and cross section values at the 
beam momentum 1628 MeV/c. 





PP 


pil7T + 


pp7r u 


d7T + 


events 


3442 


3014 


696 


80 


a mb 


21.2±0.7 


17.6±0.6 


4.48±0.20 


0.47±0.05 



3 Experimental results and discussion 
3.1 Elastic scattering 

The measured differential cross section of the elastic pp- 
scattering in c.m.s. of the reaction is shown in Fig. [1] 
The elastic cross section value given in Table [1] was calcu- 
lated as AttAq where Aq is the coefficient in front of the 
Legender polynomial Po(z) in fitting the angular distri- 
bution by the Legender expansion in the angular range 
-0.95< cos9 <0.95. 
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Fig. 1. Elastic differential cross section. The curve is result of 
the Legender polynomial fit of the -0.95< cos8 <0.95 angular 
range. 

This range was chosen in order to take into account 
a loss of events due to forward elastic scattering when a 
slow proton had a short recoil path and could not be seen 
in the bubble chamber (at the momentum less than 80 
McV/c) or was missed under the scanning. This interval 
was determined by examining the stability of the Legen- 
der coefficients with decreasing fitted range of the angular 
distribution. Our value of the total elastic cross section is 
less than that given in works [UE] approximately by 13%. 
We do not know the reason for this discrepancy: it could 
arise due to i) an insufficient control of the loss of events 
with small scattering angles or ii) an uncorrect calcula- 
tion of the millibarn-equivalent in the chamber. To con- 
trol the loss of events in the forward scattering we use the 
short-cut range for fitting of the differential cross section 
which should repair this shortcoming. Of course we have 
about 7% unmeasured events which were apportioned to 
the fraction of the fitted hypotheses for the cross section 
calculations. But it is difficult to believe that the fraction 
of missing elastic events is much larger. The missing events 
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have superimposed second tracks: this topology is not the 
elastic scattering one, except for events with disposition 
of the scattering plane strongly in the normal direction to 
the frame of a film. As concerns the millibarn-equivalent 
calculation, we would like to pay attention that the val- 
ues for the single pion production cross sections given in 
Table Q] are in a fairly good agreement with the above 
mentioned work [3]. 



3.2 The single pion production reactions. Comparison 
with OPE model 

The OPE model [TJ[2] describes the single pion produc- 
tion reactions by four pole diagrams with the ir° or ir + 
exchanges. The main evidence for pole diagram contribu- 
tions would be an observation of a peak in the momen- 
tum transfer distribution from the target particle to the 
secondary proton in the pp — > ppn process or, for exam- 
ple, to the secondary neutron in the pp —¥ pmr + process. 
Since there is no difference between final protons in the 
pp — » ppir reaction, it is difficult to separate the contri- 
bution from a certain diagram experimentally. 

Figures [2j [3] show the square momentum transfer A = 
— (p t —pf) 2 distribution for the pp — > pmr + and pp — > ppir 
reactions, where pt and pj are four-momenta of the target 
proton and one of the final nuclcons correspondingly. The 
OPE model calculations, normalized to the total number 
of experimental events, are shown there as dashed lines 
and phase space distribution as dotted ones. 
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Fig. 2. Four-momentum transfer A 2 distribution for the pp — > 
pmr + reaction. The dashed curve is the OPE calculation and 
dotted one shows shape of the phase space. 



One can see that the OPE model describes qualita- 
tively well the A 2 distribution for both reactions studied. 
It is remarkable, because only the P33 wave is taken into 
account in the intermediate irN scattering. It could be 
that this distribution is mainly sensitive to the pole dia- 
gram propagator and a more complicated structure of the 
amplitude manifests itself in other distributions. 

Figs. 2] and [5] present angular distributions of the final 
particles in the c.m.s. of the reaction as well as particle 
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Fig. 3. Four-momentum transfer A 2 distribution for the pp — > 
ppir reaction. The dashed curve is the OPE calculation and 
dotted one shows shape of the phase space. 
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Fig. 4. C.m.s. angular distributions and laboratory momen- 
tum spectra of the final particles of the pp — > ppit reaction. 
The dashed curves are predictions from the OPE model, the 
dotted curves represent the phase space. 



momentum distributions in the laboratory frame for the 
pp — > ppir and pp — > pmr + reactions correspondingly. 
One can see again that the agreement of the OPE cal- 
culations with the experimental data is fairly good. The 
qualitative agreement is also observed for other spectra, 
except for angle distributions of the final particles in the 
hclicity system. 

As noted earlier in [3], although the OPE model pro- 
vides a qualitative description of most differential spec- 
tra, it disagrees with the total cross section values for the 
pp — > ppir and pp — > pmr + reactions. It means that tak- 
ing into account the P33 (4(1232) isobar) intermediate 
state only is not enough for an adequate description of 
these reactions and a comprehensive partial wave analysis 
is needed. 
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Fig. 5. C.m.s. angular distributions and laboratory momen- 
tum spectra of the final particles of the pp — > pnir + reaction. 
The dashed curves are predictions from the OPE model, the 
dotted curves represent the phase space. 



4 Partial wave analysis 

To extract contributions from different partial waves we 
apply an event-by-event partial wave analysis (PWA) ba- 
sed on the maximum likelihood method. 

For the production of three particles with the 4-mo- 
menta qi from two particles colliding with 4-momenta k\ 
and fc 2 , the cross section is given by: 

dv= 4[ fc|^ d-P 3 (P, qi ,q 2 ,q 3 ) , P = h + k 2 , (6) 

where A is the reaction amplitude, k is the 3-momcntum of 
the initial particle calculated in the c.m.s. of the reaction, 
s = P 2 = (ki + k 2 ) 2 and d<P 3 is the invariant three-particle 
phase volume. 

The total amplitude can be written as a sum of partial 
wave amplitudes as follows: 

A = J2^)Q%... fiJ (SLJ)A 2b (i,S 2 L 2 J 2 )(s i )x 

(X 

QiZ.^S 2 L 2 J 2 S'L'J) . (7) 

Here S, L, J are spin, orbital momentum and total an- 
gular momentum of the pp system, S 2 ,L 2 ,J 2 are spin, 
orbital momentum and total angular momentum of the 
two-particle system in the final state and S',L' are spin 
and orbital momentum between two-particle system and 



the third particle with momentum q^. The invariant mass 
of two-body system can be calculated as s; = (P — qi) 2 ■ 
The multiindex a denotes all possible combinations of the 
S, i, J, L 2l J 2 , S\ L' and i, Af r (s) is the transition am- 
plitude and A 2 t,(i, S 2 L 2 J 2 )(si) describes rescattering pro- 
cesses in the final two-particle channel (e.g. the production 
of Z\(1232)). In this spin-orbital momentum decomposi- 
tion we follow the formalism given in [?l|5l[S]. The exact 
form of the operators for the initial states Q™ „ J (SLJ) 
and final states Q^" 1 „ S 2 L 2 J 2 S' L' J) can be found in 
[9] . Following this decomposition, we use spectroscopic no- 
tation 2S+1 Lj for the description of the initial state, the 
system of two final particles and the system "spectator 
and two-particle final state" . For the initial pp system, the 
states with total momenta J < 2 and angular momenta 
L = 0, 1, 2, 3 between two protons are taken into account. 
For the final three-particle system, we restrict ourselves 
in the fitting procedure by angular momenta L 2 = 0,1,2 
and V = 0,1,2. Due to nonresonant nature of the pp 
system in the energy region investigated here, there is no 
factorization of initial and final vertices, and the transi- 
tion amplitude depends on all quantum numbers which 
characterize a partial wave (index a). Moreover, due to 
contribution of the triangle singularities, the production 
parameters can be complex-valued. The best description 
was obtained with the parameterization: 



a? + a%yfs 



(8) 



where af s are real values. The a % parameter corresponds 
to a pole situated in the region of left-hand side singu- 
larities of the partial wave amplitudes. It is introduced to 
suppress the growth of the amplitudes at large s. 

Wc have also used another, more complicated param- 
etcrizations of the transition amplitude. However, we ob- 
tained either worse description of the data or similar one, 
with larger number of fitting parameters. In the latter 
case, these results serve us to determine systematical er- 
rors for the various contributions to the cross sections. For 
the ttN system in the intermediate state, we introduce two 
resonances, Z\(1232)P 33 and Roper 7V(1440)Pn. For the 
Z\(1232), we use relativistic Breit-Wigner formula with 
mass and width taken from PDG. The Roper state was 
parameterized in agreement with Breit-Wigner couplings 
found in the analysis [TU|. Let us note that the present 
analysis is not sensitive to the exact parameterization of 
the Roper resonance: only the low energy tail of this state 
can influence the data. 

For the description of the final pp interaction we use a 
modified scattering-length approximation formula: 



A 2 b(si) 



1 - \tH % 4p + *qc4 P q 2L /F(<l, rf>,L) ' 



(9) 



where multiindex /3 denotes possible combinations of kine- 
matical channel i and quantum numbers S 2 , L 2 and J 2 ; 
a^ p is the pp-scattering length and r° is the effective range 
of the pp system. The F(q, r, L) is the Blatt-Weisskopf 
form factor (it is equal to 1 for L = and the explicit form 
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for other partial waves can be found, for example, in [5]) 
and q is the relative momentum in the final two-nucleon 
system. For the S-waves this formula corresponds exactly 
to the scattering length approximation suggested in [TTJ 
[T2"] . The pp 1 Sq scattering length was fixed on value de- 
termined in the analysis [13] • In the analysis of the pp — > 
pnir + data the pn scattering length and effective range 
were fixed for S-waves at a( 1 S'o) = —23.7 fm, r( 1 S'o) = 2.8 
fm and a( 3 Si) = 5.3 fm, r( 3 5i) = 1.8 fm. 



4.1 The PWA results and discussion 

We minimized the log-likelihood value fitting the present 
data on the pp — > ppir and pp — > pnir + reactions taken 
at the proton momentum 1628 MeV/c together with ob- 
tained earlier data on pp — > ppn° [3"1I5"1[T3"] . The data [2][5] 
were taken at PNPI and measured at nine energies cov- 
ering the energy interval from 600 up to 1000 MeV. The 
high statistics data at the momentum 950 MeV/c taken 
by the Tubingen group |13j were included to fix the low 
energy region. 

The experimental data (points with error bars) and the 
results of the partial wave analysis (histograms) for the 
momentum 1628 MeV/c are shown in Figs. [5] and [71 The 
first row shows angular distributions of the final particles 
in the rest frame of the reaction and the second row shows 
effective two-body mass spectra. It is seen that our partial 
wave analysis describes these distributions rather well. 

The quality of the partial wave analysis is also demon- 
strated in angular distributions in the hclicity (third row) 
and Godfrey- Jackson (fourth row) frames. These frames 
are the rest frames of two final particle systems. In the he- 
licity frame the angle is calculated between one of the con- 
stituent particles and the spectator particle. This frame is 
mostly suitable for the investigation of cascade processes, 
when two colliding particles form a system (e.g. resonance) 
which decays into a final two-body system (e.g. another 
resonance) and a spectator. In the Godfrey-Jackson frame 
the angle is calculated between one of the constituent par- 
ticles and the beam. This system is mostly suitable to 
study production of the two-particle system due to the 
i-channel exchange mechanism. 

The initial 1 Sq partial wave provides only a small con- 
tribution to the pp —> ppir reaction at the incident proton 
momentum 1628 MeV/c. The largest contributions come 
from two P-wave initial states: 3 P 2 and 3 Pi (see Tabled]). 
The 3 Pq initial state contributes about 10% to the total 
cross section and we found notable contributions from the 
x D<i and 3 Fi partial waves. The 3 Fi partial wave inter- 
feres rather strongly with the strongest 3 Pi wave and its 
contribution is defined with a rather large error. 

All initial partial waves decay dominantly into the 
Z\(1232)p intermediate state. The contribution of channels 
with A(1232) production is varied for the different par- 
tial waves from 65 up to 90%. The strongest non-resonant 
contribution is observed from the 3 P2 initial state: here 
the transition 3 P 2 — > ( 3 P2) PP n contributed in different fits 
from 20 to 35% (from the contribution of the 3 P2 partial 
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Fig. 6. Angular distributions (a,b), effective-mass spectra of 
final particles (c,d) in the reaction rest frame and angular dis- 
tributions of final particles in the helicity (e,f) and Godfrey- 
Jackson (g,h) frames for the pp — > ppir reaction taken at the 
proton momentum 1628 MeV/c. The histograms show the re- 
sult of our partial wave analysis. 



wave). This instability appears due to a notable interfer- 
ence in this channel between the ( 3 P2) PP 'K and Z\(1232)p 
intermediate states. 

In the pp — > pnn + reaction at 1682 MeV/c, the con- 
tribution of the Z\(1232) production to the cross section 
is even stronger than in the case of the neutral pion pro- 
duction. The non-resonant pn partial waves with isospin 1 
contribute much less to the total cross section. For exam- 
ple, the transition 3 P2 ( 3 P2) pn ^ was found to be less 
than 10% from the contribution of the 3 P2 initial state. 
However we observed a notable transition 1 Sq — > ( 3 5 l i) pn 7r 
which is a dominant one for the 1 Sq initial state and an 
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Table 2. Contributions of the main partial waves to the single 
pion production reactions at 1628 MeC/c 
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Fig. 7. Angular distributions (a,b,c), effective-mass spectra of 
final particles (d,e,f) in the reaction rest frame and angular dis- 
tributions of final particles in the helicity (g,h,i) and Godfrey- 
Jackson (j,k,l) frames for the pp — > pnir + reaction taken at 
the proton momentum 1628 MeV/c. The histograms show the 
result of our partial wave analysis. 
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P, MeV 

Fig. 8. The energy dependence of the total cross section for 
the pp — > ppn° reaction. The experimental points are taken 
from articles [3ll5lll3j . The solid curve is the result of PWA; 
the long dashed curve shows the contribution from the 3 P2 
wave; the point curve from 3 Pi; the point-dashed curve from 
3 Po and short dashed curve from the 1 D 2 wave. 



appreciable contribution from the 3 P\ — > ( 1 Pi) pn 7r transi- 
tion. 

Our partial wave analysis defines relative contributions 
of the isovector waves to the total cross section of the sin- 
gle pion production processes at the energy interval 400- 
1000 MeV. Fig. [5] shows the experimental behavior of the 
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Fig. 9. Isoscalar total cross section. The dashed band rep- 
resents results of work [14]. The value calculated from eg. (1 101) 
using interpolation for o(np — > ppn^) cross section from [151 
[16] is shown as the black square with error bar. 



pp — > ppir cross section together with the result of the 
partial wave analysis and contributions of the dominant 
partial waves. It should be noted that although we use the 
data on the pp — > pnir + reaction at 1628 MeV/c only the 
found partial waves predict the total cross section at lower 
energies in a good agreement with values given in [3]). 

The present measurement of the cross section for the 
pp — > ppir® reaction together with measurements of the 
cross section for the pn — > ppix~ reaction allows us to 
obtain the isoscalar inelastic cross section by: 

a{I = 0) = 3{2<r(np — > ppir~) — <j{pp — > ppir )}. (10) 

Figure [5] shows the result of such calculation at 1628 
MeV/c using the value for the a(np — > ppir~) cross section 
interpolated from the experimental data [T31ITB] (black 
square) . The result of the work [T3] where the partial wave 
analysis of the earlier pp — > ppir data was performed to- 
gether with the np — > ppir~ measurements taken with con- 
tinues neutron beam is shown by the band. The isoscalar 
cross section calculated using the result of this analysis 
and new pp — > ppw° data at 1628 MeV is fully consistent 
with this band and is not shown here. 



5 Conclusions 

A detailed study of the differential cross section on the 
PP ~^ ppir and pp — > pni: + reactions has been performed 
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at the incident proton momentum 1628 McV/c. The shape 
of the most distributions is described qualitatively well 
by the OPE model, although it fails to describe simulta- 
neously the total cross sections for the pp — > ppir and 
pp — > pnir + reactions. 

The partial wave analysis of the single pion production 
reactions indeed reveals a dominant contribution from the 
Z\(1232)p intermediate state which explains a success of 
the OPE model. But, in addition, it allows us to obtain a 
combined description of all analyzed reactions and extract 
contributions from the transition amplitudes. 
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